Multiferroic Domain Dynamics in Strained Strontium Titanate* 
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Multiferroicity can be induced in strontium titanate by applying biaxial strain, resulting in the 
coexistence of both ferroelectric and antiferrodistortive domains. The magnitude and sign of the 
strain imposed on the lattice by design can be used to tune the phase transitions and interactions 
between these two phenomena. Using optical second harmonic generation, we report a transition 
from centrosymmetric 4/mmm phase to ferroelectric mm2, followed by an antiferrodistortive tran- 
sition to a coupled ferroelastic-ferroelectric mm2 phase in a strontium titanate thin film strained in 
biaxial tension by 0.94%. The results agree well with theoretical first principles and phase-field pre- 
dictions. Direct imaging of domains arising from the ferroelectric phase transition, and its switching 
under electric fields is demonstrated using piezoelectric force microscopy. Nonlinear optics combined 
with phase-field modeling is used to show that the dominant multiferroic domain switching mech- 
anism is through coupled 90° ferroelectric-ferroelastic domain wall motion. More broadly, these 
studies of coexisting ferroelectric (polar) and antiferrodistortive rotation (axial) phenomena could 
have relevance to multiferroics with coexisting ferroelectric (polar) and magnetic (axial) phenomena. 

PACS numbers: 77.84.Dy, 77.80.Dj, 77.80.Fm, 42.70.Mp 
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Multiferroic materials with multiple order parameters 
such as polarization, magnetization, and spontaneous 
strain lead to the coexistence of two or more of the pri- 
mary ferroic properties. As a consequence, multiferroics 
are attracting significant interest due to the possibility of 
a rich array of coupled phenomena such as ferroelastic- 
electric-magnetic, piezo-electric-magnetic, and electro- 
magneto-optic effects 0,0HSI!lil3ll. At first S lance > 
strontium titanate would appear an unlikely candidate 
for a multifcroic. Bulk SrTi03 is a cubic (m3m) cen- 
trosymmetric perovskite at room temperature; it under- 
goes a nonpolar antiferrodistortive (AFD) phase transfor- 
mation to a tetragonal point group 4/mmm at — 105 K, 
and exhibits indications of a frustrated ferroelectric (FE) 
transition at —20 K, which it never completes 0. First 
pri ncip les calculations |lfj| and thermodynamic analy- 
sis [n|, [13 have suggested however, that external strain 
can induce ferroelectricity. This prediction has been ex- 
perimentally confirmed recently in a SrTiC>3 thin film 
strained in biaxial tension [lflj . A number of fundamen- 
tal issues, however, remain to be addressed. No direct 
imaging of ferroelectric domains or its switching dynam- 
ics has been reported in this material. The multiferroic 
nature of strained strontium titanate, i.e., the coexistence 
of ferroelectric and ferroelastic domains has not received 
much attention, though it is unique in many respects, 
ft is induced by external strain. Further, unlike other 
ferroelectric-ferroelastics such as BaTiC>3 and PbTiC>3 
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which have a primary order parameter in polarization 
and secondary order parameter in strain, two indepen- 
dent primary order parameters exist: a polar ferroelec- 
tric order parameter p, and an axial antiferrodistortive 
rotation order parameter q 0, 0, 0, 0. Hence, de- 
pending on the magnitude and sign (compressive or ten- 
sile) of the applied biaxial strain, the ferroelectric transi- 
tion temperature can precede or succeed the ferroelastic 
transition temperature |l2l | . No first principles theory or 
experimental determination of the structure and polar- 
ization direction has been reported for the multiferroic 
phase. From a symmetry perspective, the SrTi03 mul- 
tiferroic domain structure has striking similarities with, 
and hence of broader relevance to other multiferroics with 
coexisting polar and axial phenomena, such as materials 
that exhibit coexisting ferroelectricity (polar property) 
and antiferromagnetism (axial property) [3|, |y, |l7J . 

This paper focuses on the domain structure and phase 
transitions in a 500 A thick SrTiOs thin film grown on 
a (110) DyScC>3 substrate with a 0.94% uniform in-plane 
biaxial tensile strain in the growth plane . Thermody- 
namic calculations predicts a ferroelectric (FE) transition 
temperature range of 277 ± 100 K, and an antiferrodis- 
tortive (AFD) transition at -111-150 K [HGJ]. First 
principles calculations were performed for this strain 
value with an Amm2 space group symmetry in the FE 
phase, and an Ima2 (#46) symmetry in the multifer- 
roic (FE+AFD) phase. In the multiferroic phase, the 
optimized ground state structure has a ferroelectric po- 
larization of p —0.18 C/m 2 in the (110) p direction in 
the growth plane, and a corresponding antiferrodistortive 
rotation of the oxygen octahedra by 5.9° about the vec- 



2 



(a)< 




8 


a 




(b)o 




o* c 


o 

o" o 






» c 


« 






o* o 


o c 


o 

o' c 




o" o 












o 

o" o 


/ 


o" C 













J 




o" c 

o 




o" c 

a 


a 

a' c 

a 




o' o 



• (c) 


V « 


® 

< ° 




V » 




(d)o 




V « 


O 

< ° 
s 


8 

p 8 


i v * 


Q* 

— <r^ 


a 


< 





® 


< c 


A 9 

a 




• < o 


\l" ® 


/ * 

< ° 


V 8 


f 




FD° 


« 


P. 8 


< o 


°V ® 


8 

< ° 


o 

V » 


< ° 




.q* 

<\ 




Dt/ c 

SE • 

+ 8 


o 




+ Sr • 


Ti o 





[010 




[1 00]p 


® 


LU 
Ll_ 





FIG. 1: (Color Online) Schematic showing Sr, O and Ti 
atomic positions and order parameters p and q in various 
phases of strained SrTiOa (a) Paraelectric i/mmm phase, (b) 
Ferroelectric mm2 phase, (c) Ferroelectric (FE) and Antifer- 
rodistortive (AFD) mm2 phase. Circles with crosses (dashes) 
are oxygen ions moving into the plane of the paper (in the 
dash direction), (d) Domains in FE+AFD phase showing an 
AFD and coupled FE+AFD domain wall. 

tor q(110) p , where the subscript p refers to the pseu- 
docubic Miller index, showing that the combination of 
the two distortions is lower in energy than either distor- 
tion alone. In this phase, there are 8 types of coexist- 
ing FE+AFD domains (see Fig. ^d)), which we clas- 
sify into 4 types of domain classes: x+,x— ,y+, and 
y— , each comprised of two types of domains. Denot- 
ing the directions x = [110] p ,j/ = [110] p , and z 
00 I ,.. the order parameters in the 4 types of domain 
classes are x+ : (p x , 0, 0, g x , 0, 0) and (^,0,0,-^,0,0), 
x- : {-p x ,0, 0,-^,0,0) and (-^,0,0,^,0,0) , 
y+ : (p y ,0,0,q y ,0,0) and (p y , 0, 0, -q y , 0, 0), and y- : 
(— Pj/,0, 0,q y , 0,0) and (— p y , 0, 0, — q y , 0, 0). These do- 
mains are separated by 180° AFD (FE) domain walls 
across which p (q) remains the same, but q (p) rotates 
by 180°; and 90° (180°) FE+AFD domain walls across 
which both p and q rotate by 90° (180°). In the pure 
FE phase, the order parameter q disappears, resulting 
in 4 types of domains, and 2 types (90° and 180°) of 
ferroelectric domain walls. 

We employ optical second harmonic generation (SHG) 
to study domain structures and phase transitions. SHG 
involves the conversion of light (electric field E") at a 
frequency uj into an optical signal at a frequency 2w by 
a nonlinear medium, through the creation of a nonlin- 
ear polarization P? u oc dij^E^ E£ . The light source for 
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FIG. 2: (Color Online) SHG polar plots at T =123 K of 
(a) if and (b) (c) if {9 = 0°) vs. T where the dark 

(hollow) squares are cooling (heating) cycle. The circles (di- 
amonds) represent Kz, y (K^l). 

SHG experiment was a ThSapphire laser with a pulse 
width of 140 femtoseconds, repetition rate of 1 kHz and 
a wavelength of 800 nm. Fundamental wave of 800 nm 
was incident from the substrate side in normal geometry. 
Its polarization direction, at an angle 9 from the y-axis, 
was rotated continuously within the plane of the film (see 
insets in Fig. EI a) and^b)). The intensity, if 3 , of the 
output SHG signal at 400 nm wavelength from the film 
was detected along either j — x,y polarization directions 
by a photomultiplier tube. The resulting polar plots of 
the SHG intensity at 123 K are shown in Fig. El A ref- 
erence study of a bare (110) DyScOa substrate without 
any film yielded no SHG intensity within the detection 
limits. Using SHG, the point group of this film at 223 
K (ferroelectric phase) was determined to be mm2, and 
the ferroelectric polarization p was found to be along the 
(110)p axes in the film plane ^2]; consistent with first 
principles calculations. A theoretical analysis of this film 
syst em y ields the following expression for the SHG inten- 
sity |U 

if = Kij sin 2 26» + K 2 ,j (sin 2 6 + K 3<j cos 2 Of 

+Ki,j (sin 2 9 + K 3 j cos 2 6) sin 20, (1) 
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where Ki j are constants given by: 

\Ks,v) \K 2tX K 2tV J U«J ' 

U^KxJ ^ - (^J • (4) 

Two of these are ratios of nonlinear coefficients (c^i / ^33 
(including its sign) and ^15/^33 (absolute value only)) 
that are independent of domain microstructure. The 
third quantity, SA x /SA y , is a microstructural quantity, 
where SA X = A x+ — A x - and 5A y = A y+ — A y - are 
the net area fractional biases in the x and y directions, 
respectively. The term SA X {8 Ay) in the observed SHG 
intensity reflect a cross-cancellation effect of signals from 
the x + (y+) domains and x — (y— ) domains resulting 
from a phase difference of it in their relative SHG fields. 
The inverse relation between K^^ x and i^3, a , stipulated 
by Eq. |2 can be seen to hold at all temperatures stud- 
ied (Fig-EIc)), showing that the symmetry is mm2 and 
the polarization is oriented along (110) p directions in the 
80-298 K range studied. These properties also show an 
anomaly in the range of T max ^250-273 K, consistent 
with the peak in the dielectric response e reported in 
Ref.|l3j|. Non-centrosymmetry and polarization, P s , ex- 
ist well up to at least Tb ^323-328 K, where finally the 
SHG signal is below our detection limit. The SHG sig- 
nal in Fig. 1(c) also shows a reproducible anomaly at 
T AFD ~ 148 K (heating) to T AFD ~ 173 K (cooling). 
This corresponds to the antiferrodistortive phase tran- 
sition in SrTiOs, consistent with thermodynamic phase 
field predictions H . The slight drop in the SHG in- 
tensity at the AFD transition could arise from a more 
effective phase cancellation of the SHG intensity due to 
doubling of the domain variants from 4 to 8 in the mul- 
tiferroic phase. It can also arise from an anomaly in the 
magnitude of across the phase transition, though the 
measured ratios of nonlinear coefficients in Fig. 2 do not 
show an anomaly. The symmetry is observed to remain 
mm2 across the AFD transition, as predicted. 

The difference in the temperatures at which the peak in 
the dielectric constant and the peak in the ^33/^31 ratio 
occurs (T max ~250 K) versus where noncentrosymmctry 
was lost (Tb ^323 K) can arise from ferroelectric relaxor 
behavior as recently reported |2(|, as well as inhomoge- 
neously strained regions. The Tb ^323-328 K transition 
may then be more appropriately termed as the Burns 
temperature [2]]]. Between T max and Tb, polarization 
in a relaxor ferroelectric is normally expected to exist as 
glassy nanopolar regions without well-defined domains or 



domain walls. This is consistent with piezoelectric force 
microscopy (PFM) at room temperature (298 K), which 
reveals that the majority of the film is featureless with 
low piezoelectric signals. Some localized regions (esti- 
mated to be in ~5-10% of the film area) however reveal 
ferroelectric domains at room temperature as seen in Fig. 
3(b), indicating inhomogeneous strains. Upon heating, 
the piezoelectric response starts to disappear partially at 
40 °C (Fig. Etc)) and completely at 60 °C (Fig.OJd)), and 
the domain features reappear on cooling (Fig. 3(e)) back 
to room temperature. The domains are partly bounded 
by ^5 nm high surface ridge-like features seen in topog- 
raphy (see Fig. 3(h)), within which they reappear after 
cooling (Fig. 3(e)). These ridge features have cracks in 
their center [22| seen typically only in SrTiOs/DyScOs 




FIG. 3: (Color Online) Topography (a), and lateral PFM 
images (b-g) at (b) 298 K (c) 313 K, (d) 333 K, (e) cooled back 
to 298 K. (f) under field of +1 kV/mm, and (g) -1 kV/mm 
in the [110] p . (h) is same as (b) but with overlapped partial 
topography from (a). 
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films thicker than 350 A. Although these films have the 
narrowest dielectric constant versus temperature peaks 
reported for SrTi0 3 or (Ba,Sr)Ti0 3 films 20], the widths 
are still about three times broader than (Ba,Sr)TiC>3 sin- 
gle crystals [2^, hence local inhomogeneous strains and 
domains can arise. As further evidence, similar films that 
are 350 A or less in thickness do not show the topo- 
graphic ridges, SHG signals or domain like features at 
room temperature. 

The bright and dark regions of PFM signal in Fig.^b) 
can be attributed to the j/± and x± domains, respec- 
tively. The domain walls (not related to topography; 
see Fig.EJh)) separating these bright and dark domains 
are parallel to the (100) p directions, which correspond to 
90° domain walls. These walls can be influenced by ex- 
ternal electric fields. Figure [jje) shows that cooling has 
produced predominantly y± domains. When an electric 
field of +1 kV/mm is applied along the +x direction, they 
transform to predominantly x± domains (Fig. 134)), and 
then partially reverse back to y± domains after applying 
-1 kV/mm along the x direction (Fig. EJg))- These re- 
sults indicate 90-degree ferroelectric domain wall motion. 
This is further confirmed by in situ SHG under in-plane 
electric fields of ±1 kV across a 1.3 mm electrode gap 
in the ±x[110] p direction, which indicates (Figures 4(a) 
and 4(b)) that both &A\ and 6A y change. This indicates 
that the electric field along the ±x direction switches do- 
mains along the ±y direction as well, indicating motion 
of 90° domain walls that couple x± domains to the y± 
domains. This mechanism is operative at all tempera- 
tures including the multiferroic phase, where SHG indi- 
cates the motion of coupled 90° ferroelectric-ferroelastic 
domain walls. 

Phase-field simulations 0, 13 were performed to 
understand the detailed mechanism of domain wall mo- 
tion under an applied electric field along +x([110] p ) in 
the multiferroic phase, as shown in Fig. Elc)-(h). One 
can clearly see that the area fraction of x+ domains in- 
creases, while that of a;— domains decreases through a 
motion of coupled 90° ferroelectric-ferroelastic domain 
walls. The pure 180° AFD walls (see Fig. Did)) within a 
single polarization domain state are degenerate in energy 
even under an electric field and remain with statistically 
equal populations (Fig. 0{f-h)). The pure AFD domain 
walls within the polarization domain states not favored 
by the external field disappear, however, primarily as the 
90° degree walls sweep through that domain state. 

In conclusion, a unique form of multiferroicity, 
namely an induced ferroelectric-antiferrodistortive phase 
in strained SrTiOa is reported. A ferroelectric tran- 
sition temperature T max ~250 K, Burns temperature 
Tb ~323 K, and an antifcrrodistortive transition at 
Tafd ~160 K have been determined using optical sec- 
ond harmonic generation (SHG). For the first time, ferro- 
electric domains were directly imaged at room tempera- 
ture and above in this material using PFM. In situ SHG 
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FIG. 4: (Color Online) Domain area fractions (J|"(0°) oc 
SAl; ly(0°) oc 8A 2 y ) as a function of applied voltage along the 
+a;([110]p) direction at (a) T=123 K and (b) T=223 K. Phase- 
field modeling of the evolution of FE ((c)-(e)) and AFD ((f)- 
(h)) domain morphologies with electric field along a;([110] p ). 

under external fields show evidence for the movement 
of coupled ferroelectric-ferroelastic 90° domain walls in 
the multiferroic phase under external fields, in agree- 
ment with phase-field modeling. From a symmetry per- 
spective, the domain structure and dynamics in this 
ferroelectric-antiferrodistortive system could have rele- 
vance to expected domain structures in a ferroelectric- 
antiferromagnetic system with two phase transitions. 
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